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An overview of the computational aspects of nonunique factorization invariants 


E A. Garcia-Sanchez 


Abstract. We give an overview of the existing algorithms to compute nonunique factorization invariants 
in hnitely generated monoids. 


1. Introduction 

In this manuscript we give a general overview of the existing procedures to deal with nonunique factor¬ 
ization invariants. These methods have gained importance since they provide batteries of examples that 
can be used to understand how to prove theoretical results (or disprove ideas that we initially thought 
would hold). The algorithms are fed from the theory and in many cases from advances in integer linear 
programming. Thus in a sense, this is a wheel: theory produces algorithms that can be used to test new 
ideas, and these yield new results. 

A semigroup is a set with a binary associative operation. If a semigroup has an identity element (an 
element that operated with any other, in both sides, keeps the element unchanged), then we say that the 
semigroup is a monoid. Let (M, •) be a monoid. An element m e M is a unit if there exists m’ e M such 
that m ■ m' - e = m' ■ m, where e is the identity element of M. A monoid is reduced if the only unit is 
the identity element. We are concerned with factorizations up to units, so we can at the very beginning 
remove the units from our monoid and consider that it is reduced. 

A monoid M is commutative ifm-m' = m'-m for all m, m' e M. All monoids in this paper are commu¬ 
tative, and thus we will adopt additive notation in the following, and will use 0 to the identity element. 

A monoid M is cancellative if whenever m + m' - m + m" for some m, m', m" g M, we have m' = m". If 
[R, -I--) is a domain, then the underlying monoid [R, •) is commutative and cancellative. As with commu¬ 
tativity, we will also assume that our monoids are cancellative. 

Thus in what follows a monoid M is meant to be commutative, cancellative and reduced. We denote 
M* =M\{0}. 

An element m in M* is said to be an atom or irreducible if whenever m = m' + m" for some m',m" e M, 
then either m' = 0 or m" = 0 (recall that we are assuming that M is reduced). Let sd (M) denote the set of 
atoms of M. We say that M is atomic if every element me M can be expressed as a sum of finitely many 
atoms. 

For a given set X, be the free monoid on X, that is, the expressions of the form T.xex AxX with 

Ax e N (N denotes the set of nonnegative integers), and all but finitely many Ax are zero. For M an atomic 
monoid, denote by Z(M) = (M)). There is a natural monoid epimorphism 

(p:Z{M]^ M, (p( Y. ^aa]= Y ^aO-. 

ass^iM) a^s^iM) 

Observe that many expressions of the form T.aEs^{M) ^aO- may correspond to the same element in M. For 
me M,we define Z(m) = (p~^im]. Every element in Z(m) is a factorization of m. For Nq M,we will write 
ZlAl) ^UmeivZ(m). 

It may happen that the cardinality of Z(m) is one for all m (and consequently qr is an isomorphism and 
M is a free monoid); in this case M is said to be a factorial monoid. It also may happen that there are 
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finitely many factorizations for every element in the monoid M, and then we say that M is a FF-monoid. 
The length of a factorization T,a£s^(M) is 'LaEs^im ^a- If for every element me M, all the lengths of its 
factorizations coincide, then we say that M is a half-factorial monoid; and if the set of possible lengths of 
factorizations are finite for every element, the the monoid is a BF-monoid (see (25l for more details and 
properties of these monoids). 

Observe that for computational aspects it is desirable that M can be described in a “finite” way, and this 
happens in the case M is an atomic monoid with finitely many atoms. In this setting, if the cardinality of 
sd (M) is e, we can identify Z(M) with As we are assuming M is cancellative and reduced, this implies, 
that any two factorizations are incomparable with respect to the usual partial ordering in Dickson’s 
lemma implies that Z(m) will have finitely many elements for any me M. 

Transfer homomorphisms allow to study the arithmetical invariants (such as sets of lengths and cate¬ 
nary degree) of Krull and weakly Krull monoids in associated auxiliary monoids. In many cases these 
auxiliary monoids are finitely generated (see (25]). So, in these cases we will have FF-monoids, and we 
will be able to determine some properties using a computer. 

Notice also that if we are assuming that M is finitely generated, then according to (351 Proposition 3.1], 
we can assume that M “lives” in x x ■ ■ ■ x Z^^. If A = {m\mg} is the set of atoms of M, then 
M = <A> = nimi \ ui,..., nf e N}. For me M the set of factorizations of m corresponds with the set 
of nonnegative integer solutions of the system of equations 

(mi I I me)[xi...Xe)^ = m, 

where the mfs are written in columns, and the last r equations are in congruences modulo di,...,ir¬ 
respectively. In order to deal with these equations in congruences we can introduce auxiliary variables 
and then project to the original ones (see for instance (35l Chapter 7]). The software Normaliz (O) can 
handle this kind of systems of equations. 

By removing equations in congruences, we then have a monoid that is, torsion free, that is whenever 
km = km’ for k a positive integer and m, m’ e M, we have m- m'. Every finitely generated commutative, 
cancellative, reduced and torsion free monoid is isomorphic to a submonoid of for some positive 
integer k (this is known in the literature as Grillet’s Theorem, see for instance |35l Theorem 3.11]). A 
monoid with all these conditions is called an affine semigroup. The set of atoms of an affine semigroup 
M is M* \ [M* M*), and it is the unique minimal generating system of M. So here minimal generators 

correspond with atoms (irreducibles). 

We will give the arithmetic invariants in the scope of affine semigroups. This does not mean that 
the some of the methods reviewed can be used in a more general scope (even in an noncomputatonal 
framework), see for instance (7l l30llMll32l . 

If z and z! are two factorizations oimeM, then the pair (z, z'] is in the kernel of the monoid morphism 
cp, which in the setting of affine semigroups with atoms {mi,..., mg} can be written as 

(fi -.N” ^ M, (p{ni,...,ng) = n\m\-\ -h ngmg. 

A presentation a of M is a generating system of ker cp = {(x, y) e N® x N® | (p[x) = (piy)}, that is, ker ^ is the 
minimal congruence containing a. 

Remark 1. Notice that from the definition of presentation, if u is a presentation for M and z, z! are two 
factorizations of me M, then there exists a chain of factorizations Zi,..., z,- of m such that 

• Zi = Z,Zr = Z', 

• for every 1 e {1,..., r -1} there exists hi, Ci e N® such that [Zi,Zi+i) = («;■ -i- c;, hi + cf with either 
{ai,bi) ea or [hi, ad e a. 

This idea actually catches the fact that ker^ is the least congruence containing a, or in other words, it is 
the reflexive-symmetric-transitive closure of a compatible with addition. 

Hence knowing a presentation of M (a generating set of ker^) allows us to know how to move from 
z to z', and consequently it will be a fundamental tool in the study factorizations of elements in affine 
semigroups. This is the case of tame and a»-primality. 
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Recently it has been shown that some invariants are related to the calculation of the set of factoriza¬ 
tions of a principal ideal, and if the monoid is full, there are specific procedures that significantly speed 
the process. 

For numerical semigroups there are particularizations of the procedures based mainly in Apery sets, 
which avoid the use of linear integer programming, and work well for small generators. We will describe 
them when applicable. 

This manuscript is meant to give a state of art of the implementations existing for the calculation of 
nonunique factorization invariants. We will simply explain the theory that supports these procedures, 
but will not describe deeply the functions used. We have implemented everything that is described here 
in the GAP (HU) package numericalsgps {(IU)> and thus it is part of this package (see the manual of the 
package for a description of the functions, examples and mode of operation). The reader interested in a 
full description and implementation of the algorithms can have a look at the source code available either 
in the GAP web page, or for the developmentversion in https: //bitbucket. org/gap-system/numericalsgps 
(the files containing the functions described here for numerical semigroups are in catenary-tame .gi 
and contributions .gi; those for affine semigroups are in affine.gi, both in the folder gap). The 
package tests availability of other packages ( (2411211 [HI l2Zl 121) that interact with 4ti2 ([T]), Normaliz 
(lU) and Singular (HU). Depending on this availability, the package will use an specific method for the 
calculations. So we had to implement in some cases up to four methods for the same invariant depending 
on the extra software used (this is why there are several files with prefix af f ine-extra in the gap folder). 

2. Presentations 

Redei proved in (SU that every finitely generated commutative monoid is finitely presented. In our 
setting this means that every affine semigroup admits a presentation with finitely many elements. Since 
then, many alternative and shorter proves have been published. We recall here one of these approaches. 

Let t be an unknown and k be a field. For M an affine semigroup, define the semigroup ring k[M] = 
where addition is performed component-wise and multiplication follows the rule = 

^m+ml 

Assume that {mi.mg} is a generating system of M. Herzog in (2U proves that u is a presentation of 

M if and only if the ideal Im = (^“ - I ici, b] e a], where Im is the kernel of the ring homomorphism 
induced by 

k [xi.Xf] —" k [M], X; '->• f . 

Observe that for n = (ui,..., nj;), we can write t'” as and in this way we can see k[M] as a subring 

of k [ q,..., ffc]. In particular, we can compute a presentation of M by using elimination: we start with the 
ideal (xi-t'”',...,Xe-r'”®) £ k[xi,...,Xe, q,..., tfc], and then eliminate the variables q,..., to obtain/ m- 

Example 1. Let us compute a presentation of M= ((2,0), (0,2), (1,1), (2,1)) with singular, HU . 

> ring r=0,(x,y,z,t,u,v),Ip; 

> ideal i = (x-u~2,y-v~2,z-u*v,t-u*v~2); 

> eliminate(i,u*v); 

_[l]=yz2-t2 

_ [2]=xt2-z4 
_ [3]=xy-z2 

This means that Im = [yz^ - - z^, xy - z^), and in light of Herzog’s correspondence, the set 

{((0,1,2,0), (0,0,0,2)), ((1,0,0,2), (0,0,4,0)), ((1,1,0,0), (0,0,2,0))} 
is a presentation for M. 

A minimal presentation of M is a presentation that cannot be refined to another presentation of M, that 
is, it is minimal with respect to set inclusion (it turns out that it is also minimal with respect to cardinality; 
see (SU Corollary 9.5]). 
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Example 2. The presentation in Example[T]is not minimal. If we want to obtain a minimal presentation 
with singular additional work is needed. 

> ring r=0,(x,y,z,t,u,v),(wp(2,2,2,3),lp(2)); 

// ** redefining r ** 

> ideal i = (x-u**2,y-v**2,z-u*v,t-u*v**2); 

> ideal j=eliininate(i,u*v); 

> minbase(j); 

_ [1]=xy-z2 

_ [2]=yz2-t2 

Given m e M, we define as the graph with vertices Z[m) and zz' is an edge if z-z' 5^ 0 (dot product). 
An element m is a Betti element of M if the graph is not connected. We will denote by Betti{M) the set 
of Betti elements of M. 

The sets of vertices of the connected components of are also known as ^-classes of Z(m). The 
following method can be used to produce all minimal presentations (up to arrangement of the pairs and 
symmetry) of M; see for instance |35l Chapter 9]. 

• For all me M, if is connected, then set Om - t>- If not, let i?i,..., be the different 5?-classes 

ofZ(m). Consider any tree r with vertices For each i e {1,..., r} take riCRi. Set a fn = 

{[Zi,Zj) I RiRj is an edge of T} (for instance, one might take am - {izi,Z 2 ), (zi, Z3), ..., (zi, Zr)}). 

• The set a = [JmEMt^m is a minimal presentation of M. 

It follows that the set of Betti elements of M has finite cardinality. And that the cardinality of a (any) 
minimal presentation is LfoeBetti(M) (ncc(V/,)-!), where ncc(V jj) stands for the number of connected com¬ 
ponents of Vj,. This formula holds for every atomic monoid having the ascending chain on principal 
ideals ((71 Corollary 1]). 

Example 3. Fet M be as in Fxample[I] Since any presentation contains a minimal presentation, we have 
that Betti(Af) £ {(2,4), (2,2), (4,4)}. We use the GAP ((H]) package numericalsgps (Hg)) to calculate the 
5?-classes of each of these elements. 
gap> RClassesOfSetOfFactorizations( 

FactorizationsVectorWRTList([4,4],[[2,0],[0,2],[1,1],[1,2]])); 

[ [ [ 0, 0, 4, 0 ], [1, 0, 0, 2 ], [ 1, 1, 2, 0 ], [ 2, 2, 0, 0 ] ] ] 
gap> RClassesOfSetOfFactorizations( 

FactorizationsVectorWRTList([2,4], [[2,0],[0,2], [1,1],[1,2]])); 

[ [ [ 0, 1, 2, 0 ], [ 1, 2, 0, 0 ] ], [ [ 0, 0, 0, 2 ] ] ] 
gap> RClassesOfSetOfFactorizations( 

FactorizationsVectorWRTList([2,2], [[2,0],[0,2], [1,1],[1,2]])); 

[ [ [ 1, 1, 0, 0 ] ], [ [ 0, 0, 2, 0 ] ] ] 

It follows that Betti(AT) = {(2,2), (2,4)} (this also follows from Fxample|2). 

The function FactorizationsVectorWRTList either uses (El, or if available (5| or [T1 through the 
packages Normalizinterface ((27|) or either 4ti2gap ((24]) or 4ti2Interf ace ((261). 

We will see that the catenary degree and the Delta sets are “ruled” by a minimal presentation. 

3. Apery sets 

Let M be an affine semigroup generated by {mi ,..., mg}. Let me M. The Apery set of m in M is the set 

Ap(M, m) = [m' e M\m' -m^ M}. 

Apery sets can be defined in a more general setting. If our monoid fulfills the ascending chain condition 
on principal ideals, then every for every m' e M there exits unique {w,k] e Ap(M, m) x N such that m' = 
km + w (see 13). 
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If M is a numerical semigroup, then the cardinality of Ap(M, m) has exactly m elements. Moreover, if 
b e Betti(M), then b = mi + w with i e{2,...,e} and iv £ Ap(M, mi) \ {0} (see for instance EEl Proposition 
8.19]). 

As minimal presentations are crucial for studying factorizations, the above paragraph implies that 
Apery sets are also important in our study particularized to the numerical semigroup setting. 

4. Graver bases 

Let M he an affine semigroup, generated by {mi,..., me). 

We have seen that a minimal presentation is a minimal generating system of kerip as a congruence. It 
turns out that ker^ is not only a congruence, but an affine semigroup itself, and thus it admits a unique 
minimal generating system, which we denote by J^(M). It follows easily that J^{M) corresponds with the 
pairs {x,y) = ({xi,...,Xe),{yi,...,ye)) e x N'^\{{0,0)} that are minimal (with respect to the usual product 
order) solutions of 

(mi I • • • I mel - mi I • • • I - me) (x I y) ^ ^ 0, 

because if (x, y) e ker^, thenximiH-tXgme = yimiH-tyemg. Moreover, there exists (xi,yi),..., (Xs,y 5 ) e 

J^{M) such that (x, y) = (xi, yi) H- 1 - (x^, y^). That is, every pair of factorizations of the same element can 

be expressed as a sum of factorizations of some specific elements. Indeed, we will say that m e M is 
pnmfime if there exists x,yeZ(_m) such that (x, y) 

In particular, J^{M) is a presentation of M, though in general with a lot of redundancy. 

Notice that if e, is the ith row of the identity ex e matrix, then (e/,e/) e J^(M) for all i e (1,..., e}. 

On define the order (xi,...,Xe) E (yi,...,ye) if for all i e {l,...,e}, x/y, > 0 and |x, | < |y, |. Also, for 
X e set x"*" and x~ to be the unique elements in such that x = x'*' - x“ and x"*" • x“ = 0. It turns out 
that X E y if and only if (x"*", x“) < (y''', y~) (usual partial ordering). 

Let iT be a subgroup of Z®. A Graver basis of iZ is a set of minimal nonzero elements of H with respect 
to E. 

Notice that the set of integer solutions of 

(mi I ••• I me)x^ = 0 

defines a subgroup Hm of Z®. In fact (x, y) e ker^ if and only if x - y e Hm (this is a rephrasing of the 
necessity condition in (35j Proposition 1.4]). From a Graver basis G of Hm we can easily compute 

J^iM) = {(x+,x“) I xeG} u{(ei,e,) j i e |l,...,e}}. 

Example 4. Let us go back to M in Examples[I]and[2j 

gap> GraverBasis4ti2(["mat".TransposedMat([[2,0], [0,2],[1,1],[1,2]])]); 

[ [ 1, 0, -4, 2 ], [0, 1, 2, -2 ], [ 1, 1, -2, 0 ], [1, 2, 0, -2 ] ] 

The output of 4ti2 does not print an element and its negation. Hence a Graver basis of Hm consists in 8 
elements and J^(M) has 8+4 elements. 

We will see that some nonunique factorization invariants depend on the factorizations of the primitive 
elements of M. 


5. Block monoids 

Let G be an Abelian group. And let gi,..., gfc ^ G- A zero-sum sequence is an expression of the form 
uigi + —h nicgk = 0. The length of this sequence is ni + —h n^. We say that a zero sum sequence is min- 
imalii there is no other zero-sum sequence n'^gi + ■•• + n'^gk - 0 such that b ^ {n[,...,n'^) < (ni,...,nG- 
The set of zero-sum sequences is clearly a monoid, actually it can be identified as a submonoid of and 
it is generated by the minimal zero-sum sequences (indeed it is a full affine semigroup). We will denote 
the set of zero-sum sequences in gi,..., g^ by ^ ({gi,..., g^}). 


6 


P. A. garcia-sAnchez 


Since G is an Abelian group, it is then isomorphic to x • • • x x t} for some di,...,dr,l eN. Hence 
we can identify the elements with elements in x x Z^^ xZK Hence corre¬ 

sponds with the set of nonnegative integer solutions of the system of r -i- / equations and k unknowns 

(gl I ■ • • I gk)x ^ 0 £ Zrfj X • • • X Zd^ X z' 

(the first r equations are in congruences modulo di,...,dr, respectively). The set of solutions of this 
system of equations can be computed via Normaliz ((5)). 

The Davenport constant is the supremum (in this setting maximum) of the lengths of minimal zero- 
sum sequences. 

Examples. We can compute the block monoid associated to Z^ in the following way using numericalsgps. 
gap> m2:=[[0,l],[1,0],[1,1]];; 

gap> a:=AffineSemigroupC"equations",[TransposedMat(m2),[2,2]]);; 
gap> GeneratorsQfAffineSemigroup(a); 

[ [ 0 , 0 , 2 ], [ 0 , 2 , 0 ], [ 1 , 1 , 1 ], [ 2 , 0 , 0 ] ] 

Observe that we are omitting (0,0) and that the second argument of Af f ineSemigroup is a matrix whose 
columns are the elements in (Z^)* and a list indicating the equations that are congruences with the re¬ 
spective modules. The Davenport constant in this case is 3. 

Many factorization properties of monoids can be derived from the factorization properties (or bounded 
in some cases) of the block monoid of their class groups (see (25]). This is why these affine semigroups 
are relevant in the study of nonunique factorization invariants. 


6. Denumerant and maximal denumerant 

We have already mentioned that for an affine semigroup M and me M, the set Z(m) has finitely many 
elements. The denumerant of m is precisely the cardinality of Z{m). There is a wide amount of literature 
devoted to the study of denumerants of elements in numerical semigroups, indeed few formulas are 
known, and just for some particular families of monoids ((33) is a nice reference for the reader interested 
in this topic). 

Of course the bigger an integer in a numerical semigroup is, the larger is its denumerant, and thus it 
is not bounded. What is indeed astonishing is that the maximal denumerant is bounded for numerical 
semigroups. The maximal denumerant of m in M is the number of elements in Z(m) with maximal 
length. If M is a numerical semigroup, set the maximal denumerant of M as the maximum of the ma xi mal 
denumerants of elements of M. Bryant and Hamblin give in (6] a procedure to compute the maximal 
denumerant of any numerical semigroup. 

Example 6. The semigroup <3,5,7) has maximal denumerant 2. 

gap> s:=NumericalSemigroup(3,5,7);; 

gap> MaximalDenumerantOfNumericalSemigroup(s); 

2 

gap> List(IntersectionC[0..100],s), 

> x->Length(FactorizationsElenientWRTNumericalSemigroup(x,s))); 


[ 1, 

1, 1, 

1, 

1, 

1, 1, 

2, 

1, 

2, 2, 

2, 

3, 

2, 3, 

3, 

3, 

4, 4, 

4, 

4, 5, 5, 

5, 

6, 

6, 

6, 7, 

7, 

7, 

8, 8, 

9, 

9, 

9, 10 

10, 11, 11 

, 12, 12, 12 

, 14, 13, 14 

:, 15, 

15, 

16, 

16, 

17, 

17, 

18, 

19, 

19, 

20, 

20, 

21, 

22, 

22, 

23, 

24, 

24, 25, 

26, 

26, 

27, 

28, 

29, 

29, 

30, 

31, 

31, 

33, 

33, 

34, 

35, 

35, 

37, 

37, 

38, 

39, 40, 

41, 

41, 

43, 

43, 

44, 

46, 

46, 

47, 

48, 

49, 

50, 

51, 

52, 

53, 

54, 

55 ] 
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7. Length based invariants 

Let M be an affine semigroup generated by {mi,..., me). Take me M and x = (xi,..., Xg) e Z{m). Recall 
that the length of x is defined as 

|x| = Xi + ••• + Xg. 

The set of lengths of factorizations of m is 

L(m) = {|x| I X e Z(m)}. 

Since Z(m) has finitely many elements, so has L(m). This means that affine semigroups are BF-monoids. 
Recall that a monoid is half factorial if the cardinality of L(s) is one for all s e S. This concept was 
introduced for domains in [41] . 

From Remarkdjit easily follows that M is half factorial if and only if for every [a, b) in a minimal pre¬ 
sentation of M we have \a\ = \b\ (see (381)- Thus we can determine whether or not an affine semigroup is 
half-factorial. 

Example 7. In Example[2j since ((1,2,0,0), (0,0,0,2)) belongs to a minimal presentation of M, we deduce 
that M is not half-factorial. 


7.1. Elasticity. One of the first nonunique factorization invariants that appeared in the literature was the 
elasticity (introduced in [40]). It was meant to measure how far is a monoid from being half factorial. 
Take m in an affine semigroup M. The elasticity of m, p[m), is defined as 

^ ^ supL(m) 
p{m)= . I . 

mmL(m) 

Since L(m) has finitely many elements, the supremum in the numerator is indeed a maximum. The 
elasticity of M is defined as 

p(M) = sup{p(m) I m e M]. 

It is not hard to show (see [38]) that 


f \a\ 

p{M)^max\ — 
[ \b\ 


{a,b)eJ{M) 


Hence hy computing a Graver basis of Hm we can calculate the elasticity of M. Philipp in his thesis, 
and published later in [3l] , did a great improvement in the computation of the elasticity: he showed 
that we only have to consider elements [a, b) e J?(M) with a h and with minimal support (indices of 
nonzero coordinates). These elements are known in the literature as circuits, and we can use QT] Lemma 
8.8] to calculate them by means of determinants. The calculation of a Graver basis is in general a hard 
problem, while computing determinants is affordable. Thus Philipp theoretical contribution became a 
considerable speed up in combination with Eisendbud and Sturmfels method for computing circuits. 


p(M) = max 


\b\ 


[a, b) circuit of ker^ 




Example 8. Let us compute p 

gap> in: = [[0,0,l] , [0,1,0] , [0,1,1] , [1,0,0] , [1,0,1] , [1,1,0] , [1,1,1]] ;; 
gap> a:=AffineSemigroupC'equations",[TransposedMat(m),[2,2,2]]);; 
gap> ElasticityOfAffineSemigroup(a); 

2 


7.2. Delta sets. Another way to measure how far we are from half factoriality, is to determine how distant 
are the different lengths of factorizations. This is the motivation for the following definition. 

Let as above m be an element in the affine semigroup M. Assume that L(m) = (Zi < • • • < Ir). Define the 
Delta set of m as 


A(s) — {I2 — h,...,lr ~ lr-l}> 
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and if r = 1, A(m) = 0. The Delta set of M is defined as 

MM]= U 

meM 

So, the higger A(Af) is, the farther is M from begin half factorial. 

Recall that (x, y) e ker(p if and only if x - y £ Hm- Indeed, it is not hard to show that Hm is generated as 
a group by the differences of the pairs in a presentation of M. From this, one can prove that 

minA{M) = gcdA(AD 

( (25l Proposition 1.4.4]). 

By using the idea expressed in Remark[T] it can be shown that the maximum of the distances between 
lengths of factorizations is reached in a Betti element of M ([TT] Theorem 2.5]): 

maxA(M) = max{maxA(fi) | b e Betti(AD}. 

This gives us an interval where the elements in A(M) must be, but it is far from being a procedure to 
compute the whole set A(M). 

For numerical semigroups, it is known that the sets of distances between consecutive lengths of factor¬ 
izations are eventually periodic (US]) and a bound for this periodicity is given. This bound was improved 
in [20]. Hence we can compute the Delta sets of the elements up to this bound (a dynamic version of this 
procedure is resented in O). The problem is that this bound can be huge. 
gap> s:=NumericalSemigroup(701,902,1041); 

<Numerical semigroup with 3 generators> 
gap> DeltaSetOfNumericalSemigroup(s); 

[ 1, 2, 3, 4, 5, 6, 11, 17 ] 

gap> DeltaSetPeriodicityBoundForNumericalSemigroup(s); 

313436 

Recently in (22] a procedure that runs as fast as Euclid’s extended algorithm has been presented for nu¬ 
merical semigroups with embedding dimension three (and not symmetric, though the algorithm seems 
to work also for symmetric numerical semigroups). 

O’Neil in (29| gives new tools based on Hilbert functions that probably will yield procedures for the 
computation of A(M) for an arbitrary affine semigroup M. 

8. Distance based invariants 

Observe that length based invariants cannot describe the behavior of factorizations in half-factorial 
monoids. To measure how spread are the factorizations, we first need a distance. 

For X = (xi,..., Xg), y = (yi,..., ye) e N®, define the infimum of x and y as 

xAy= (min{xi,yi},...,min{Xp,yp}) 

(if we think in additive notation and x and y are factorizations of an element, then x A y translates to 
greatest common divisor). 

The distance between x and y is defined as 

d(x,y) = max{|x-(xAy)|,|y-(xAy)|} 

(equivalently d(x, y) = max{|x|, |y|| - |xAy|). 

8.1. Catenary degree. We start with an example that illustrates the idea of catenary degree. 

Example 9. The factorizations of 66 e (6,9,11> are 

Z(66) = {(0,0,6), (1,3,3), (2,6,0), (4,1,3), (5,4,0), (8,2,0), (11,0,0)}. 

The distance between (11,0,0) and (0,0,6) is 11. 
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In the above picture the factorizations are depicted in the top of a post, and they are linked by a “catenary” 
labeled with the distance between two consecutive sticks. On the bottom we have drawn the factoriza¬ 
tions removing the common part with the one on the left and that of the right, respectively. So we have 
linked (11,0,0) and (0,0,6) with a chain of factorizations, and every two consecutive nodes in the chain 
are at most at distance 4. This is in fact the best we can do in this example. We are not caring about the 
length of the sequence, but about how closer are two consecutive elements in the chain. 

Let M be an affine semigroup, and take me M. Let x, y e Z(m) and let AT be a nonnegative integer. An 
Al-cham joining x and y is a sequence xi,...,X/t e Z(m) such that 

• Xi = X, Xfc = y, 

• for all I e {1,..., fc- 1}, d(x;,X;+i) < N. 

The catenary degree of m, denoted dm], is the least N such that for any two factorizations x, y e Z(m), 
there is an Al-chain joining them. The catenary degree of M, c{M), is defined as 

c(M) = sup{c(m) I m e M}. 

The calculation of dm) can be performed in the following way. We consider the complete graph with 
vertices the factorizations of m, and edges labeled with the distances between their ends. Then we pick 
an edge with the largest label, and if it is not a bridge, then we remove it. We keep doing so, until we arrive 
to a bridge. The label of this bridge is c(m). 

Example 10. As an illustration of the above procedure, consider 77 e S = <10,11,23,35). In the following 
figure we see that we can remove the edge with label 6, meaning that in order to go from (0,7,0,0) to 
(2,1,2,0) we can first go to (2,2,0,1) and then to (2,1,2,0), and the distances in this walk between two 
consecutive nodes are less than 6. Then we remove the edge labeled with 5. But we cannot remove the 
edge joining (1,4,1,0) and (0,7,0,0) since it is abridge (we can remove the other labeled with 3). 



Thus the catenary degree of 77 is 3. 




10 


P. A. garcia-sAnchez 


Observe that in Remark [TJ we obtained chains joining any two factorizations of the same element, 
just using translations of elements in a presentation. Since distances are not translation-sensitive, we 
only have to care on how to go from the first component to the second in a relation in a presentation. It 
follows (see HOl) that 

c(M) = max{c(i>) | b e Betti(M)}. 

This gives a computational procedure to compute the catenary degree of any affine semigroup M. 

Example 11. Let us recover Example^ M = ((2,0), (0,2), (1,1), (2,1)). We already know that Betti(AD = 
{(2,2), (2,4)}. 

gap> a:=AffineSemigroupC[2,0],[0,2],[1,1],[1,2]);; 
gap> gens:=GeneratorsQfAffineSemigroup(a); 

[[0,2], [1,1], [1,2], [2,0]] 
gap> betti:=BettiEleinents0fAffineSemigroup(a); 

[[2,2], [2,4]] 

gap> List(betti,b->FactorizationsVectorWRTList(b,gens)); 

[ [ [ 1 , 0 , 0 , 1 ], [ 0 , 2 , 0 , 0 ] ], 

[ [ 2 , 0 , 0 , 1 ], [ 1 , 2 , 0 , 0 ], [ 0 , 0 , 2 , 0 ] ] ] 

gap> List(last,CatenaryDegreeOfSetOfFactorizations); 

[ 2, 3 ] 

gap> CatenaryDegreeOfAffineSemigroup(a); 

3 


So far we do not know of a procedure to compute the (finite) set {c(m) | m e M}. It is known that for 
numerical semigroups, the catenary degree is also eventually periodic, but unfortunately no bounds for 
this periodicity are known ([8]). For half-factorial monoids it can be shown (see |23l Theorem 2.3]) that 

{c(m) I m G M} = {c(m) | m e Betti(AT)}. 

For numerical semigroups, in light of Section[3](see also O Corollary 3]), 

c(M) = max{c(m) | me {m 2 ,...,me}-)-(Ap(M, mi)\|0})}, 

and so in this setting it is not needed to compute Betti(M). 

8.2. Monotone, equal and homogeneous catenary degrees. We can obtain different flavored catenary 
degrees if we impose conditions on the definition of AT-chain. For instance if we enforce the chain of 
factorizations to have nondecreasing lengths we obtain the definition of monotone catenary degree. 

We can also ask the lengths to be all equal, and then we have equal catenary degree. 

Finally we can also impose that the lengths in the chain are not larger than the maximum of the lengths 
of the ends of the chain, obtaining in this way the homogeneous catenary degree. 

Let M be an affine semigroup, M Q N^. Let me M,we write (m, 1) e for the element with the first 
coordinates the coordinates of m and last coordinate equal to 1 (we have appended a 1 at the “end” of 
m). Assume that M is minimally generated hy (mi,..., mg}. Set 

AT'^^ = ((mi,l),...,(m„l)> 

and 

M^°“-((mi,l),...,(me,l),(0, D). 

Both and AT^°™ are half-factorial monoids (see (23) ). It is easy to prove that the equal catenary degre 

of M corresponds with the catenary degree of AT®'?, and that the homogeneous catenary degree of M is 
precisely at^'’"'. This provides a way to compute both homogeneous and equal catenary degrees. 

In order to compute the monotone catenary degre of AT, it can be derived from (31] that we have to look 
at the projections in the first k coordinates of the primitive elements of AT^'’™ (see ES Chapter 3]), and 
then take the maximum of the monotone catenary degrees of these elements. The monotone catenary 
degree of m is the maximum of the equal and adjacent catenary degree of m, where the adjacent catenary 
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degree of m is defined as follows: let L(m) = {Zi < • • • < 1^}, and for every i e denote by Zi. (m) the 

set of factorizations of m with length the adjacent catenary degree of m is the maximum of the distances 
d(Z;,.,Z;;^j), Z e r- 1}. 

Example 12. Let us use numericalsgps to compute the catenary degrees of (10,17,24,31,43). 
gap> s:=NumericalSemigroup(10,17,24,31,43); 

<Numerical semigroup with 5 generators> 
gap> MinimalGeneratingSystem(s); 

[ 10, 17, 24, 31, 43 ] 

gap> CatenaryDegreeQfNumericalSemigroup(s); 

6 

gap> HomogeneousCatenaryDegreeOfNumericalSemigroup(s); 

11 

gap> MonotoneCatenaryDegreeOfNumericalSemigroup(s); 

11 

gap> EqualCatenaryDegreeOfNumericalSemigroup(s); 

11 

8.3. Tame degree. Assume that M is an affine semigroup generated by {mi,..., mg}, and let min M and 
X e Z{m). If there exists ui,..., ng e N such that m - (L^=i t^i) ^ then there must be y = (yi. • • • > ye) ^ 
Z(m) such that y - (ni,..., Hg) e N®. We want to know the smallest possible distance at which we can find 
such an y. This is the idea of tame degree. We are mostly interested in the case L;=i m, = mj for some 
je{l,...,e}. 

The tame degree of m with respect to mi, t[m, mi), is the least nonnegative integer t such that for every 
z e Z(m), there exists z! c Z(m) with z! - e; e (or in other words, the Zth coordinate of z! is nonzero) 
and d(z, z') < t. The tame degree of M with respect to m;, t(M, m;), is the supremum (maximum in this 
setting, [To]) of all the tame degrees of the elements of m; + M with respect to m;. 

The tame degree of M, t(M), is the maximum of the tame degrees of S with respect to all the atoms 
(affine semigroups are tame and locally tame, (23). The tame degree of M can be computed by means 
of the tame degrees of the primitive elements of M ([101). Recently, a faster approach has been described 
in [21]. Set Jii = Minimals<Z(mi + M) and M, = {(piz) \ z e Mi\. By Dickson’s lemma, JZi and Mi have 
finitely many elements. Moreover, 

t(M, mi) = max{t(m, m,) | m e M,}. 

In [37] there is a procedure to compute Mi (indeed the set of expressions of any ideal of M, not just 
principal ideals). By using [5] or [T] (or any integer linear programming package) we can also compute 
this directly as in the following example. 


Example 13. Let us compute the set Mi for M = ((2,0), (0,2), (1,1), (1,2)>. We need to find the expressions 
in (2,0) + M. This corresponds with the (x, y, z, t) e such that 






fxl 

'2 

0 

1 

1] 

y 

0 

2 

1 

2] 

z 

J J 


= (2,0) + 


'2 0 1 li 
0 2 12 


y' 

z' 

[t'j 


for some (x', y', z!, t') e N^. And this is a system of two equations and eight unknowns. We use the package 
4ti2gap to solve this. 

gap> m:=[[2,0,l,l,-2,0,-l,-l],[0,2,l,2,0,-2,-l,-2]]; 

[ [ 2, 0, 1, 1, -2, 0, -1, -1 ], [ 0, 2, 1, 2, 0, -2, -1, -2 ] ] 
gap> problem: = ["mat" ,in, "rhs", [[2,0]] ,"sign", [ [1,1,1,1,1,1,1,1] ]] ; 

[ "mat", [ [ 2, 0, 1, 1, -2, 0, -1, -1 ], [ 0, 2, 1, 2, 0, -2, -1, -2 ] ], 

"rhs", [ [ 2, 0 ] ], "sign", [ [ 1, 1, 1, 1, 1, 1, 1, 1 ] ] ] 









12 


P. A. garcia-sAnchez 


gap> ZSolve4ti2(problem); 


rec ( 

zhom 

; = 

[ [ 

1, 

0, 

0, 2, 

0, 

0, 

4, < 

3 ] 

, [ 

0, 

0, 

4, 

0, 

1 

[ 

1, 

2, 

0, 

0, 

0, 

0, 

0, 

2 ], 

[ 0, 

0, 

1, 

0, 

0, 

0, 

1, 

0 ] 

3 


[ 

0, 

0, 

0, 

1, 

0, 

0, 

0, 

1 ], 

[ 0, 

1, 

2, 

0, 

0, 

0, 

0, 

2 ] 

3 


[ 

0, 

0, 

2, 

0, 

1, 

1, 

0, 

0 ], 

[ 0, 

1, 

0, 

0, 

0, 

1, 

0, 

0 ] 

3 


[ 

0, 

0, 

0, 

2, 

1, 

2, 

0, 

0 ], 

[ 1, 

0, 

0, 

0, 

1, 

0, 

0, 

0 ] 

3 


[ 

1, 

1, 

0, 

0, 

0, 

0, 

2, 

0 ], 

[ 0, 

0, 

0, 

2, 

0, 

1, 

2, 

0 ] 

] 


zinhom 

: = 

[ [ 

0, 

0, 

4, 

0, 0, 

0, 

0, 

2 ] 

, [ 

0, 

0, 

2, 

0, 

0, 

1 

[ 

1, 

0, 

0, 

0, 

0, 

0, 

0, 

0 ], 

[ 0, 

0, 

0, 

2, 

0, 

2, 

0, 

0 ] 

] 

) 


This in particular means that 


Z((2, 0) + M] = {(0,0,4,0), (0,0,2,0), (1,0,0,0), (0,0,0,2)} 

+ ((1,0,0,1),(0,0,4,0),(1,2,0,0), (0,0,1,0), (0,0,2,0),(0,1,0,0), 10,0,0,2), (1,0,0,0),(1,1,0,0)> 

AndthusMinimals<Z((2,0) + M) ^ {(0,0,2,0), (1,0,0,0), (0,0,0,2)}. Hence Mi ^ {(2,2), (2,0), (2,4)}. 


If M is a full affine semigroup (for instance in the case of block monoids), then the elements in can 
he computed using (4) Corollary 3.5]. In this case Mi is the set of minimal nonnegative integer solutions 
of 

(mi I ••• I me]x^ > rrii. 

Example 14. Let us compute as explained in (21] the tame degree of SS (Z^). 

gap> c: = [ [ 0, 0, 1 ], [ 0, 1, 0 ], [ 0, 1, 1 ], [ 1, 0, 0 ] , [ 1, 0, 1 ], 
[ 1 , 1 , 0 ], [ 1 , 1 , 1 ]]; 

gap> a:=AffineSemigroupC'equations",[TransposedMat(m),[2,2,2]]);; 
gap> TameDegreeOfAffineSemigroup(a); 

4 


For numerical semigroups we have a similar behavior as in the catenary degree. The tame degree is 
reached in an element that has to do with Apery sets ((5J Theorem 16])): 


t(M) = max 


t(m) 


m G 


{mi,...,me} + 


e \] 

\jApiM, mi) \m 

H=i J 


For small generators, the above formula is faster than computing minimal factorizations in principal 
ideals (or if we do not have software to solve linear Diophantine equations over the set of nonnegative 
integers at hand). 


9. W-PRIMALITY 

Let M be an affine semigroup. Define on M the following binary relation: m <m m' if m' -me M. This 
relation is an order relation (the translation of divisibility to additive notation). We say that me Mis prime 
if whenever m <m m' + m" for some m', m” e M, either m <m m' or m <m m”. Any prime element must 
be an atom. But it may happen that no atom is prime (this holds in any nontrivial numerical semigroup). 
The to-primality is meant to determine how far is an element from being prime. 

The (x)-primality of m in M, denoted (x)[m), is the least positive integer N such that whenever m <m 
fli + ••• + for some ai,...,a„e M, then m <m ai^ + --- + ai^ for some {q.ijv] £ {!,...,«]. 

According to this definition an element is prime provided that its w-primality is one. 

Notice that by definition, m <m m' if and only if m' is in the principal ideal m + M. Hence princi¬ 
pal ideals play a fundamental role in the computation of a»-primality (as in the calculation of the tame 
degree). Indeed in (H Proposition 3.3] it is shown that 

(x){m) - max{|x| | x e Minimals<(Z(m -i- Af))}- 
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In (H] the above formula together with the algorithm presented in (37) is used to compute the oi-primality 
of an element in an affine semigroup. One can also proceed as in Example [T3] and use for instance 
Normaliz or 4ti2. 

The omega primality of M, if M is minimally generated by {mi,..., mg}, is defined as a){M) as the max¬ 
imum of {a)imi),...,a)ime)}. Note that ei(m) with m running in M is not upper bounded in general. 

Example 15. According to Example (131 for M = ((2,0), (0,2), (1,1), (1,2)), we have ai({2,0)) = 2. Let us 

double check it with the numericalsgps package. 

gap> a:=AffineSemigroupC[2,0],[0,2],[1,1],[1,2]);; 

gap> OmegaPrimalityOfElementlnAffineSemigroupC[2,0] ,a); 

2 

gap> OmegaPrimalityOfAffineSemigroup(a); 

4 


Eor numerical semigroups we obtain a similar construction as for the tame degree (as expected, since 
we are using roughly the same elements in the calculations). In (H Remarks 5.9] it is shown that if we 
are looking for minimal factorizations in Z(m -i- M), then we only have to search for factorizations of the 
elements of the form m+w with weAp[M, rrii) for some i e {1,..., e}. In (3] an improved method that also 
uses Apery sets is given (its actually the one that uses the package numericalsgps; see contributions .gi 
in the package). 

Example 16. Let us compare the timings for S = (10,17,24,31,43). 
gap> s:=NumericalSemigroup(10,17,24,31,43);; 
gap> QmegaPrimalityQfNumericalSemigroup(s);time; 

11 

13 

gap> a:=AsAffineSemigroup(s);; 

gap> OmegaPrimalityOfAffineSemigroup(a);time; 

11 

3654 

(The timings are in milliseconds.) 

If the generators are larger, then the principal ideal approach is better. 
gap> s:=NumericalSemigroup(201,223,357);; 
gap> OmegaPrimalityOfNumericalSemigroup(s);time; 

75 

32245 

gap> a:=AsAffineSemigroup(s);; 

gap> OmegaPrimalityOfAffineSemigroup(a);time; 

75 

1934 
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